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ABSTRACT 

As part of the Herschel Multi-tiered Extragalactic Survey we have investigated the rest-frame 
far-infrared (FIR) properties of a sample of more than 4800 Lyman Break Galaxies (LBGs) in the 
Great Observatories Origins Deep Survey North field. Most LBGs are not detected individually, 
but we do detect a sub-sample of 12 objects at 0.7<z< 1.6 and one object at z = 2.0. The ones 
detected by Herschel SPIRE have redder observed NUV-U and U-R colors than the others, while the 
undetected ones have colors consistent with average LBGs at z > 2.5. The UV-to-FIR spectral energy 
distributions of the objects detected in the rest-frame FIR are investigated using the code CIGALE to 
estimate physical parameters. We find that LBGs detected by SPIRE are high mass, luminous infrared 
galaxies. It appears that LBGs are located in a triangle-shaped region in the ^fuv vs. LogLpuy 
diagram limited by Apuv = at the bottom and by a diagonal following the temporal evolution of 



4 



BurgarcUa ct al. 



the most massive galaxies from the bottom-right to the top-left of the diagram. This upper envelop 
can be used as upper limits for the UV dust attenuation as a function of Lpuv- The limits of this 
region are well explained using a closed-box model, where the chemical evolution of galaxies produces 
metals, which in turn lead to higher dust attenuation when the galaxies age. 

Subject headings: galaxies: evolution — galaxies: formation — galaxies: high-redshift — infrared: 
galaxies — ultraviolet: galaxies 



1. INTRODUCTION 

The safest way to estimate the total star forma- 
tion rates (SFR) of galaxies is to consider the energy 
budget involving far-ultrav iolet (FUV) and far-infrared 
(FIR) measurements (e.g.. IBuat fc Xulll996f ). But, be- 
cause only a small number of individual Lyman Break 
Galaxies (LB Gs) have been detecte d in the FIR/submm 
range (e.g.. iChapmaii et"al] l2000t iChapman fc CasevI 
120091 : iSiana et al.ll2009| ). we need to observe this type of 
galax ies at lower redsh i fts to understand their dust emis- 
sion. IBurgarella et al.l (|2007D detected dropout galaxies 
at z ~ 1 at 24 /im with Spitzer. But, the dust luminosi- 
ties estimated from the rest frame 12 fim flux density is 
far from the peak of the dust emission and could provide 
biased SFR estimates. 

We observe in the FIR a sample of LBGs at 
0.7 < z < 1.6 (FUV dropouts) and at 1.6 < z < 2.8 (near- 
UV or NUV dropouts). We use the SPIRE instru 
ment (jGrifSn et all |2010[) on Herschel (iPilbratt et al 
|2010[) with observations from HerMES (jOliver et al 
[Mot lOliver et al.l l2010bf Fl. This is the first opportu- 
nity to estimate directly the dust luminosity (or upper 
limits) of unlensed LBGs from FIR data. We assume 
flm — 0.3, r^A = 0.7, and Hq = 71 kms^^ Mpc^^ and use 
AB magnitudes throughout. 

2. DATA 
2.1. LBG samples 

Observations of the Great Observatories Origins Deep 
Survey North (GOODS-N) were secured as part of the 
Herschel Multi-T iered Extragalactic Survey (HerMES, 
lOliver et al.l [2010l ) by GALEX in FUV and NUV. We 
define two samples of galaxies in two redshift ranges cor- 
responding to FUV dropouts and NUV dropouts (Fig.[T]). 

The ph otometry is performed with IRAF daophot ii 
(jStetsonl [1987) in the NUV, and in the FUV with the 
NUV coordinates. Using ADDSTAR, the completeness 
is estimated to 80% down to m[F UV1 = 24 9 and 
to[NUV] = 24.2. We use cigale (INoll et al.l 120091) 
to build models in the range < z < 3 and delimit 
the regions corresponding to LBGs at 0.7<z<1.6 and 
1.6 < z < 2.8 in the color-color diagr ams H The sample is 
cross-correlated with the i?-selected iCapak et al.l (j2004l ) 
multiwavclcngth catalog over 0.4 deg^, with a search ra- 
dius of 1", producing a catalog of 86,768 entries (46,076 
with GALEX data, 47,450 with optical data and 6,979 
objects with both) . The U -band data were collected us- 
ing the Kitt Peak National Observatory 4-m telescope. 
The V-, R-, /-, and z-band data were collected us- 
ing the Subaru 8.2-m telescope and Suprime-Cam instru- 
ment with 5(7 limiting magnitudes of [/ = 27.1, B = 26.9, 



denis . burgarellaQoamp.frJ 
scbasticn.hcinisSSoamp.tr|l 
^ http://hGrmes.sussex.ac.uk 

^ Note that the code accounts for IGM attenuation. 



V = 26.8, R = 26.6, / = 25.6 and z = 25. 4. Photometric 
redshifts arc computed using Le Phare (jArnouts et al.l 
|2002: Ilbcrt ct al. 2009) and we use them in addition 
to spe ctroscopic redshifts from iBarger. Cowie. fc Wand 
(|2008| ). In the redshift range 0.7<z<1.6, we compare 
photometric to spectroscopic redshifts (784 galaxies) and 
find CTz / (1+z) = 0.036 while we have (24 galaxies) CTz / 
(l+z) = 0.125 at 1.6 < z < 2.8. 

From the lower redshift catalog, 27 objects (4%) have 
photometric redshifts at z < 0.7, 696 objects (94%) are in 
our redshift range of 0.7 <z < 1.6 and 17 objects (2.3%) 
are at z > 1.6. For the higher redshift catalog, 223 objects 
(5%) have photometric redshifts at z < 1.6, 3859 objects 
(94%) are in our redshift range of 1.6<z<2.8, and 25 
objects (0.6%) are at z > 2.8. 

2.2. Matching with SPIRE 

We cross-correlate the LBG s amples with the HerME S 
multi- wavelength catalog from iRoseboom et al.l (j2010D . 
It uses 24 fim sources as a prior, so we restrict ourselves 
to the deep GOODS-N MIPS region of 19' x 12', i.e., 
0.063 deg^ Within this region we have 260 low-z LBGs 
and 1558 high-z. This HerMES catalog contains 1951 
24 /im sources as an input to the association process. A 
search radius of 1" is chosen, since the 3.6 /xm positions 
are good to that accuracy. Given the source density of 
the HerMES catalog this should return < 1% spurious 
matches. We find 86 matches between this list and our 
LBG samples, with 63 matches to the low-z sample and 

23 matches to the high-z sample, with ~ 2 (3.2%) and 
^ 11 (47.8%), respectively, expected by chance. 

The HerMES catalog gives SPIRE measurements for 
all MIPS sources. To define confident SPIRE detections, 
we require fiux densities higher than 7.6, 9.2 or 10.4 mJy 
at 250, 350 and 500 /xm respectively. This flux limit cor- 
responds to 2 X CTfaint, the cl ipped map con f usion noise, 
where CTfaint is estimated by iNguven et al.l (j2010[) from 
the map variance after removing pixels brighter than 
SCTconf (the raw map confusion noise). Imposing these 
more reasonable flux density limits gives 14 low-z candi- 
dates and one high-z candidate. 

We dcflne SPIRE detections (listed in Table [T]) when- 
ever tlie_SNR^^s_Jargcr than 3 a nd their "puri ty" index 
(see IRoseboom eFal . 2010; Bris bin et al.ll20iot ) is larger 
than 0.2. Note that according to this definition there 
arc no detections at 500 /xm. The "purity" of the SPIRE 
flux density is estimated from the ratio of this source's 

24 /im flux density to the 24 /xm one smoothed with the 
SPIRE beam at this position. For no "pollution", the 
purity is 1 and decreases when there are possible contri- 
butions by other sources. We additionally inspected the 
images around each candidate to check for contamina- 
tion by neighbors. Based on the purity index, we exclude 
2 objects at 0.7<z<1.6. All but one of the sources on 
the final sample (actually the highest redshift LBG) have 
SNR> 6 (see Table [1]). 
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Fig. 1.— The upper panel shows the selection for 0.7<z< 1.6 dropouts: m[FUV]-m[NUV] > 0.7(m[NUV]-S)+2.5; m[FUV]-m[NUV] > 
2.5; m[NUV] - B < 1.8; and SNR(NUV) > 3. This corresponds to 740 objects, with an average SNR{NUV) = 8.6. The lower panel shows 
the selection for 1.6 < z < 2.8 dropouts: m[NUV] - B > 1.25(B - z) + 1.5; m[NUV] - B > 1.5; m[NUV] - B < 5.0; B - z < 0.4; and 
SNR(B) > 3.0, which corresponds to 4107 objects, with a sample average SNR(i3) = 9.7. The LBGs are shown as blue dots and the 
SPIRE-detcctod LBGs by big (red) boxes. Both LBG samples arc plotted on each panel. Yellow diamonds arc stars. 



The final 0.7 < z < 1.6 HerMES LBG sample contains 
12 objects and the 1.6 < z < 2.8 sample contains 1 ob- 
ject. Although it is difficult to determine what fraction 
of z ^2 LBGs should be detected in the FIR, we can 
use the Spitzer 24 /im data to estimate how many z ~ 1 
LBGs are expected to be detected at 250 /im. From the 
present sample, the mean ratio S250/S2A ~ 54 ± 11 for 
the SPIRE-cletcctcd LBGs and so, the present 250 /im de- 
tection limit of 7.6 mJy correspond to 140 /iJy at 24 /im. 



At this level, and based on the 24 /im catalog, we would 
expect 20 z 1 LBGs to be detected with SPIRE at 
250 /im. We find 12 objects, i.e., 60% of the expectation. 
T his is consistent w i th the completeness at this flux level 
in lRoseboom et al.l (|2010( ). 



2.3. S ED fitting 
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Dust luminosities (^ir = L[8— 1000 fim] ) and ot her pa- 
rameters are estimated using CIGALE (jNoll et al.ll2009l FI. 
CIGALE performs a Bayesian analysis to estimate param- 
eter by fitting models to the UV-to-submm SEDs. One 
can select among two single stellar population libraries 
and several IR models/templates. An AGN component 
can also be added to estimate the AGN fraction (con- 
tribution by a potential AGN to Lir). The parameters 
of the dust attenuation law can be modified and CIGALE 
allows for two separate stellar populations with a multi- 
phase dust treatment. 

We use th e VLA 1.4 GHz radio data from 
iMorrison et al.l (|2010[ ) which provides data down 
to a RMS noise of 3.9 /xJy p e r bea m. The stellar 
emission is based on |Maraston' ('20 Q5I), wh ile the dust 
emission is based on [Dale & Heloul |2002) templates. 
CIGALE provides dust luminosities Ljr, while FUV 
luminosities Lfuv arc derived at Arcst ~ 153 nm and arc 
defined through the quantity vL^. An energy budget 
is performed during the fit, and the maximum value 
allowed for Ljr has to be consistent with the energy 
moved by dust grains from the UV-optical range to the 
FIR range. 

Individual SEDs with the best models selected by 
CIGALE are shown in Fig. [3] for the five LBGs with radio 
data and for the (z = 1.9) LBG. Observational and phys- 
ical parameters are given in Table [TJ Fig. [3] shows that 
CIGALE is able to fit the SEDs from the FUV to the radio 
successfully. The average FUV luminosity of the SPIRE 
detected objects is log((LFUv/Lo)} = (10.7 ± 0.2), 
and most of them are therefore UV-Luminous Galaxies. 
Their average dust luminosity is log((LiR/L0)) = (11. 9± 
0.1) and their stellar average mass is log((M^,/Mp)) = 
(11.0 ± 0.5). A comparison with iMagdis et all (|2010f ) 
shows that the average stellar mass is similar to the av- 
erage stellar mass (log((M*/M0)) = 11) of z ~ 3 LBGs 
detected at A = 8 /um. The dust temperatures Td are es- 
timated for a few objects by fitting modified black bodies 
with an emissivity index of 1.5. We find that two of the 
low-z LBGs have Td ~ 24 K, while the high-z ULIRG 
has Td = 52K (see Table [I]). The two low-z LBGs have 
quite low temperatures compared with ULIRGS/SMGs, 
but may be typical of star-forming galaxies in general, 
while the high-z one is similar to e.g. Arp220. 

3. DUST ATTENUATION OF LYMAN BREAK GALAXIES 

Our LBG SEDs in the rest-frame UV are in very 
goo d agreement with the rest- frame UV spectra of 
the lShaplev et al.l ()2003| )'s composite spectrum at z~3. 
Whatever the redshift range, LBGs seem to present 
the same starburst characteristics in the rest-frame UV, 
which is expected, given that similar rest-frame color se- 
lection criteria are used. However, the colors become 
redder when the objects are detected at 24 /im by Spitzer 
and even redder when they are detected at 250 fim by 
Herschel SPIRE. This can be interpreted as being due 
to higher d ust attenuations (s ee Fi g. [5]), and is con- 
sistent with iBurgarella et al.l (|2007f ). who found that 
both^pztzer-detected and undetected LBGs have about 
the same stellar population ages, but the latter are more 
extinguished. 

^ http://www.oamp.fr/cigale 



The AGN fraction determined by CIGALE is al- 
ways consistent with zero, except for the galaxy 
J123633.2-K620834.3, for which a 20% AGN contribution 
to the FIR dust luminosity is suggested. Most of the IR 
excess due to a potential AGN should be in the MIR, 
while the SPIRE fiux is e xpected to be dominated b y 
the starburst component (jHatziminaoglou et al.ll2010() : 
CIGALE accounts for both. 

The vast majority of the dropout galaxies are not 
detected by SPIRE with the pr esent de tection limits. 
If LB Gs were to follow IMeurer . Hcck man. fc Calzettil 
(|1999D relation (relating UV-attenuation to dust emis- 
sion), what fraction of them would be detectable? We 
can use the {U — V) color as a proxy for the rest-frame 
m[FUV] — m[NUV]. Simulations suggest that this color 
provides (3 with an accuracy better than 0.01 for power- 
laws f\ oc A*^ and /3 = —2, —1 and 0. In detail, we can 
estimate; 

/3 = log[/,(FUV)//,(NUV)]/log[AFuv/ANUv] - 2.00 
= 5.391og[5^(FUV)/5^(NUV)] - 2.00. (1) 

From the /3 values, we estimate log(LiR/I/FUv) for the 
LBGs and, after evaluating Lfuv from the filter clos- 
est to AFuv=0.15nm (as a function of the redshift), 
we determine Ljr. These values of Lir can be trans- 
formed into L250 and 5'25 n using the fol l owing calibra- 
tions computed from the iDale fc Heloul (|2002[ ) models 
with 1.0<a<2.5, i.e., star-forming galaxies with sim- 
ilar properties t o LB Gs. We have also checked that 
ICharv fc Elba3 (pOOl models are consistent with our 
calibration. By fitting polynomials as a function of red- 
shift we find: 

3 

logiiR = + log(^^-) U , (2) 

4=0 

with C = (Co,Ci,C2,C3) = (1.168, -1.166, 0.565, 
-0.091) for A = 250^lm, (1.689, -1.426, 0.582, -0.089) 
for A = 350 /xm, and (2.336, -1.605, 0.556, -0.078) for 
A = 500 /im. 

Assuming the lMeurer. Heckman. fc Galzettil (|1999() re- 
lation allows to estimate rough order of magn itude 
Lir/Lfuy for our LBGs (|Burgarella et al.l l2009t) . At 
250 /im we should have detected 10 LBGs from the sam- 
ple undetected in the mid-IR and far-IR, three LBGs 
from those detected at 24 /im, and 2 LBGs from those 
detected at 250 /im. 

For the LBGs not detected by Herschel, we can es- 
timate upper limits of log(LiR) for each SPIRE band 
and we use the lowest of the three values as the the 
final upper limit on Lir/Lfuv- From log(iiR/LFUv) 
we can estimate ^fuv- Fig. |4| suggests that the maxi- 
mum level of attenuation depends on Lfuv ~ the most 
UV-luminous LBGs yield a lower maximum ^fuv than 
less UV-luminous ones. This is true for upper lim- 
its as well as for detections. The two higher redshift 
lensed LB Gs detected in F IR or in s ub-mm (the "Cos- 
mic Eye" , ISiana et al.ll2009l and cB58. lSiana etalll200l 
also comply with this upper boundary, if we correct for 
the amplification. As do the two u nlensed LBGs MM8 
(IChapman et al.ll2000D and MMDll (jChapman fc CasevI 
|2009[ ). We stress that most of the upper limits should 
populate the area below the observational limits but none 
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TABLE 1 

Observational data and physical parameters deduced by cigale. The column origin redshift is "specz" if the redshift is 

spectroscopic or the 3(t uncertainty otherwise. 



ID z origin ;u[NUV] U Sn S250 S3B0 S1..1GH., log J/, SFR logLm ,4fuv logtpuv Ti 
redaliift [AB] [AHj [^Jy] [mjy] [mjy] [mjjj [Mq] [Mpyr-'] [Mq] M [Le] [K] 



J123624.6+620610.2 


0.75 


0.11 


22.87 


± 


0.04 


23.25 


± 


0.04 


93.2 


± 


8.3 


8.0 


± 


0.9 


















10.0 


± 


0.1 


23.5 ± 0.048 


11.1 


± 


0.039 


1.96 


± 


0.14 


10.4 


± 


0.6 




J123547.4+621005.9 


0.81 


0.24 


24.07 


± 


0.07 


25.33 


± 


0.24 


86.3 


± 


7.6 


11.3 


± 


0.9 


12 


1 


± 


1 


5 








10.2 


± 


0.1 


34.8 ± 0.030 


11.3 


± 


0.007 


3.68 


± 


0.28 


9.9 


± 


1.4 




J123724.8+620938.5 


0.81 


0.13 


22.80 


± 


0.04 


23.59 


± 


0.06 


249.0 


± 


6.7 


10.3 


± 


0.8 


















10.7 


± 


0.1 


39.4 ± 0.090 


11.3 


± 


0.094 


2.42 


± 


0.28 


10.5 


± 


0.9 




J123633.2+620834.9 


0.93 


specz 


22.41 


± 


0.03 


22.84 


± 


0.03 


779.0 


± 


7.9 


22.4 


± 


1.2 


19 


2 


± 


1 


8 


50.8 


± 


10.2 


11.1 


± 


0.2 


119.7 ± 0.040 


11.9 


± 


0.027 


2.84 


± 


0.18 


10.9 


± 


1.1 


24.3 


J123714.4+622112.3 


0.94 


specz 


24.20 


± 


0.08 


25.17 


± 


0.21 


212.0 


± 


4.6 


9.2 


± 


1.2 


12 





± 


2 










10.7 


± 


0.2 


43.9 ± 0.059 


11.4 


± 


0.050 


3.46 


± 


0.32 


10.1 


± 


1.3 




J123624.4+620836.3 


0.95 


specz 


23.22 


± 


0.05 


23.83 


± 


0.07 


695.0 


± 


11.5 


21.1 


± 


2.8 


22 


9 


± 


7 











11.3 


± 


0.2 


116.4 ± 0.182 
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Fig. 2. — In this {U — R) vs. {m[NUV] — U) diagram, blue points are LBGs which are undetected at 24 fim, green dots arc those detected 
at 24 fim and red boxes are those detected in the FIR with Herschel. We see on the marginal distributions that LBGs get redder with 
increasing maximum wavelength of detection, especially in the U-R color. This is expected if IR-bright LBGs arc more attenuated than 
IR-faint LBGs. 



would prevent us from detecting LBGs that would have 
larger FUV dust attenuations than the one suggested by 
t he present data for a given ipuv- 

iReddv et al.l (|2010( ) suggest that objects with a lower 
UV luminosity at 1.5 < z < 2.6 have lower bolometric lu- 
minosities than UV-bright galaxies which, in turn, may 
suggest lower dust attenuations if we account for the re- 
lation between bolometric luminosities a nd dust attenu - 
ations. This is in agreement with Bouwcns et al.l (|201Clf ) 
at much higher redshifts. But, when compu ting Apjjv 
from the stacked points of iReddv et al.l (120101) . we find a 
trend simila r to ours f or mos t of the points. ICarilli et al.l 
(|2008[ ) and iHo et al.l (|2010[ ) estimate the UV dust at- 
tenuation by comparing radio-based star formation rates 
to UV-based ones using a stacking analysis. We show 



Fig- in the two points corresponding to different log Lfuv 
and th ey agree with the above trend. [Burgarell a et al.l 
(|2006( ) reached a similar conclusion using Spitzer 24 ^m 
observations of a sample of LBG at z~l. Finally, we 
have divided our z 1 sample in two sub-samples as a 
function of hpjjv and stacked them in the 250/xm im- 
age (Fig. 2]). We find the same trend, again. Note, 
howeve r, that a ccounting for the error bars, our stack- 
ing and IReddv et al, (2010[ ).'s are both consistent with a 
constant Apuv- Also, whenever the selection is not fully 
complete, one may miss objects in the regions where the 
brightest sources (UV and IR) lie. This effect due to the 
inhomogene ous back ground produces holes in the stack- 
ing (see e.g. iBavouzc t 200a). For this very reason, the 
z 1 stacked points are only considered as lower limits 
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Fig. 3. — Observed spectral energy distributions of LBGs superimposed on best fit models (S^fmjy] vs. logA[/tm]). Note tfiat 
J123633.2+620834.9 very likely hosts an AGN, as suggested by the CIGALE SED fitting. 



The first LBGs individually detected in the FIR by Herschel 
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and we were not able to stack z ^ 2 LBGs. 

To understand the origin of this effect, we b uild a sim- 
ple closed box model (see iPagel et al.|[l997D . assuming- 
several exponentially decreasing star formation histories 
= *oe"*/^, with T = 0.1, 1 and lOGyr. We as- 
sume a mass of cold gas Aigas that forms stars following 
a Salpeter initial mass function, and thus produce heavy 



elements. Mgas evolves as follows: 



c;MgasM = -*(t) + S(t), 



(3) 



where E{t) is the mass ejected by stars at the end of their 
lifetime. 

The oxygen abundance Zq can be estimated as 



Zf 



-Po ln[l - a(l 



-t/r 



(4) 



where po is the oxygen yield and a is the fraction of mass 
kept in st ellar remnant s . We estimate an empirical rela- 
tion from iReddv et all ()2010[ ) that links 12 + log(0/H) 
to Lir: 

log(LiR/LFuv) = 1.67 (12 + log(0/H)) - 12.72. (5) 



IReddv et al.l (|2010f l's objects are in the rcdshift range 
1.5 <z< 2.6, so very close to ours and with metallicities 
in the range 8.2 < 12/log(0/H) < 8.8 which corresponds 
to ages in the range of a few 100 Myrs to 5 Gyrs. So, 
strictly speaking, our models are extrapolations for low 
r's but are in the good range for r = 10 Gyrs. Finally, 
from Lm we comp ute Apuy usin g the relationship from 
[Biu'garella. Buat, fc Iglesias-Para mo (2005]). 

Fig. m shows that this simple closed box model follows 
the same trend as our LBGs in the diagram. The initial 
mass of gas is set to log(Mgas/M0) = 10.5 to explain 
the low redshift LBGs and log(Mgas/MQ) = 11.0 (not 
plotted) for the high rcdshift LBGs. This is in agree- 
ment with the mass of cold gas predicted by models 
(e.g., ILacev et al.l I201H) . We find that star- formation 
timescales of r < 1 Gyr seem to be ruled out by this 
model. 

4. CONCLUSIONS 

We have selected two samples of LBGs at z ^1 and 
z~2. For the first time, we can put constraints on the 
dust emission and therefore the dust attenuation of LBGs 
directly from rest-frame FIR measurements of individual 
LBGs observed with Herschel-SPIRE. Two main conclu- 
sions can be drawn from this analysis. 

We detected 12/260 - 4.6% and 1/1558 ~ 0.06% of the 
LBGs at 0.7 < z < 1.6 and 1.6 < z < 2.8, respectively. AU 



the other LBGs arc undetected by SPIRE, and their dust 
attenuation is lower than the detected LBGs.. However, 
we have to account for the fact that the limits depend on 
ipuv- 

Secondly, the maximum dust attenuation in the FUV 
decreases as UV luminosities increase. 
Other points of interests are: 

• The dropout selections presented in this paper are 
very efficient (~^95%) at detecting galaxies in the 
redshift range 0.7 < z < 2.8. 

• The rest frame UV SEDs of the two dropout sam- 
ples are similar to higher redshift LBGs. 

• CIGALE is able to model the observed SEDs from 
the FUV to the radio and we provide the derived 
physical parameters. 

• The stellar masses of these IR-bright dropout 
galaxies are of the same order as the stellar masses 
of IR-bright LBGs observed in IRAC and MIPS 
bands. 



We propose that all LBGs lie in a triangle-shaped 
region in the ^fuv vs. LogLpuv = diagram 
limited by dust- free (small and/or young) galaxies 
to the bottom and by the locus for evolving most 
massive galaxies to the top. 



SPIRE has been developed by a consortium of insti- 
tutes led by Cardiff University (UK) and including Univ. 
Lcthbridge (Canada); NAOC (China); CEA, DAMP 
(France); IFSI, Univ. Padua (Italy); lAC (Spain); 
Stockholm Observatory (Sweden); Imperial College 
London, RAL, UCL-MSSL, UKATC, Univ. Sussex 
(UK); and Caltcch/JPL, IPAC, Univ. Colorado (USA). 
This development has been supported by national 
funding agencies: CSA (Canada); NAOC (China); 
CEA, ONES, CNRS (France); ASI (Italy); MCINN 
(Spain); Stockholm Observatory (Sweden); STFC (UK); 
and NASA (USA). The data presented in this paper 
will be released through the Herschel Database in 
Marseille HeDaM (http: / /hedam.oamp.fr/ HerMES ) 
This work mafes iise oi TOP CAT 



http://www.star.bristol.ac.uk/~mbt/topcat/ 
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